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Abstract. We observed with the camera FORSl on the 
VLT (UTl, ANTU) 50 randomly selected fields and ana- 
lyzed the cosmic shear inside circular apertures with di- 
ameter ranging from 0.5 to 5.0 arc-minutes. The images 
were obtained in optimal conditions by using the Service 
Observing proposed by ESO on the VLT, which enabled us 
to get a well-defined and homogeneous set of data. The 50 
fields cover a 0.65 square-degrees area spread over more 
than 1000 square-degrees which provides a sample ideal 
for minimizing the cosmic variance. Using the same tech- 
niques as in Van Waerbeke et al. 2000, we measured the 
cosmic shear signal and investigated the systematics of the 
VLT sample. We find a significant excess of correlations 
between galaxy ellipticities on those angular scales. The 
amplitude and the shape of the correlation as function of 
angular scale are remarkably similar to those reported so 
far. 

Using our combined VLT and CFHT data and adding the 
results published by other teams we put the first joint con- 
straints on r^o and erg using cosmic shear surveys. From 
a deduced average of the redshift of the sources the com- 
bined data are consistent with erg ~ 0.59lg Qgrio"" "*^ (for 
a CDM power spectrum and F = 0.21), in excellent agree- 
ment with the cosmological constraints obtained from the 
local cluster abundance. 



Send offprint requests to: mellier@iap.fr 

* Based on observations obtained at the Very Large Telescope 
UTl (ANTU) which is operated by the European Southern 
Observatory (program 63.O-0039A). 
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1. Introduction 

The weak lensing by large scale structures (hereafter cos- 
mic shear) probes the projected mass density along the 
line-of-sight, regardless of the distribution of light. This 
could yield powerful constraints on the evolution of the 
dark matter power spectrum, the cosmological parame- 
ters and the mass/light bias as a function of the angu- 
lar scale and the redshift (see reviews from Mellier 1999, 
Bartelmann & Schneider 2001 and references therein). 

Four teams announced recently the first measure- 
ments of a significant variance of the shear (7^) on 
scales below 30 arc-minutes (Van Waerbeke et al. 2000, 
Bacon et al. 2000a, Wittman et al. 2000, 

Kaiser et al. 2000). All the measurements are in remark- 
able agreement, and it is worth to mention that these anal- 
yses were done from independent data sets obtained on 
various telescopes and by measuring galaxy shapes with 
different tools. Moreover the signal is consistent with the 
theoretical expectations of current cosmological models, 
and does not show strong contamination by residual sys- 
tematics. Though it is not yet possible to break the degen- 
eracy between the normalization of the power spectrum, 
(Tg, and the density parameter fig from a measure of (7^) 
alone (Bernardeau et al. 1997, Jain & Seljak 1997), some 
cosmological models can already be rejected (like the 
COBE-normalized SCDM which predicts too much power 
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at small scale, at the 5-a level). These encouraging re- 
sults show that cosmic shear is close to be a mature tool 
and offers a neat complementary approach in cosmology 
compared to CMB anisotropies, SNIa, or galaxy catalogue 
statistics. 

However, more quantitative and reliable constraints 
on cosmological models are still under way. The need 
for a high confidence level on the cosmic shear sig- 
nal demands an excellent control of the systemat- 
ics and a significant reduction of the cosmic vari- 
ance (Van Wacrbckc ct al. 1999). The former issue is 
discussed at length by Van Waerbeke et al. 2000 and 
Bacon et al. 2000a on real data, and it is also addressed 
in Kaiser 2000, Erben et al. 2000 and Bacon et al. 2000b 
using simulated data. They have carried out convincing 
tests and simulations in order to demonstrate that grav- 
itational weak distortion is measurable down to slightly 
below the percent level, without critical issues regarding 
residuals from the PSF corrections. 

In order to convince that the signal is not contami- 
nated by some unexpected systematics, it is still impor- 
tant to diversify the observing conditions, and produce 
more independent data sets with different systematics. Of 
equal importance, the field- to-field variations, caused by 
source and lens clusterings and by cosmic variance, are an 
additional source of noise. The four studies mentionned 
above cover a total area of about 5 deg^ but only totalize 
28 uncorrelated fields. Clearly, one needs to increase first 
the number of uncorrelated fields, rather than the total 
solid angle of the survey, to explore the field-to-field vari- 
ation of the cosmic shear amplitude and to minimize the 
cosmic variance. 

This work explores a new data set obtained with the 
faint object spectrograph FORSl (Appenzeller et al 1998) 
mounted on the first VLT (UTl/ANTU) at Paranal (pro- 
gram 63.O-0039A; PI: Mellier). These data provide in- 
dependent measurements of the cosmic shear signal, in 
addition to the four previous studies, but with a totally 
new technology telescope and a new observing mode. 
In contrast to previous telescopes the VLT is continu- 
ously operated in an active optics mode and offers the 
possibility to observe in service mode. Thanks to its 
large collecting power, the VLT can observe many sin- 
gle fields very quickly which permits to do deep pointings 
on many different lines-of-sights. This new strategy, com- 
plementary to the one we adopted on the CFHT data 
(Van Waerbeke et al. 2000), gave us an interesting large 
complementary sample of 50 VLT/FORSl uncorrelated 
fields for cosmic shear. The analysis of these fields is dis- 
cussed below. They are then combined with other existing 
results to give the first constraints on {Qo,as) derived from 
cosmic shear. 

The paper is organized as follows: Section 2 is a de- 
scription of the data and the observations. Details on data 
reduction are given in Section 3. The analysis of the fields 
and the measurements of the shapes of galaxies are devel- 



oped in Section 4. In Section 5 we discuss the systematics. 
Discussion and combination with other results in a cos- 
mological context is given in Section 6. 

2. Data and Observations 

The targets were selected from three very large sky ar- 
eas close to the North and South galactic poles with the 

following criteria: 

1. no stars brighter than 8th magnitude inside a circle of 
1 degree; around the FORS field (to avoid light scat- 
tering); 

2. no stars brighter than 14th magnitude inside the FORS 
field; 

3. no extended bright galaxies in the field. Their luminous 
halo could contaminate the shape of galaxy located 
around them; 

4. no rejection of over-dense regions, where clusters or 
groups of galaxies could be present from visual inspec- 
tion of the Digital Sky Survey, and no preferences to- 
wards empty fields. This criterion avoids biases toward 
under-dense regions with systematically low value of 
the convergence; 

5. minimum separation between each pointing of at least 
5 degrees in order to minimize the correlation between 
fields; 

6. galactic latitudes lower than 70° in order to get enough 
stars per field for the PSF correction. 

We selected 50 FORSl fields, each covering 6.8' x 6.8'. The 
total field of view is 0.64 deg^ and the pointings are ran- 
domly distributed over more than 1000 deg^. Among those 
50 fields, two have been selected from the STIS parallel 
data sample, for cosmic shear analysis on very small scale, 
and two are in common with the WHT sample observed 
by Bacon et al. 2000a ^. 

The observations are in I-band data only and were ob- 
tained with FORSl on the VLT/UTl (ANTU) at the 
Paranal Observatory. They were carried out in service 
mode by the ESO staff during Period 63 (March 1999 
to September 1999) with the standard imaging mode of 
the spectrograph. It turns out that our program is per- 
fectly suited for service mode, both from a practical point 
of view, since they are easily set up during service mode 
periods on UTl, and from a scientific point of view, since 
we have a guarantee to get a complete set of data within 
our specifications. 

FORSl is equipped with a 2048x2048 thinned CCD, back- 
side illuminated. The pixel size is 24 /im, providing a scale 
of 0.2 arc-second in standard mode of the instrument^. Re- 
markably, no fringe pattern has been detected on FORSl 
with this filter, so I-band observations enable one to get 

^ Wc plan to check the reliability of the measurements by 
comparing the analyses done by each group on these common 
fields 

^ see http://www.eso.org/instruments/forsl/index.html 
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optimal image quality on each field. Since we only re- 
quested non-photometric nights, with seeing lower than 
0.8 arc-seconds, the schedule of our 50 fields over the 
semester was indeed quite easy. The total exposure time 
(36 minutes) was computed in order to reach / = 24.5, 
which corresponds to a galaxy number density per fields 
of about 30 gal arcmin^^. The expected average redshift 
of the lensed sources is {z) ~ 1. Each pointing was split 
into at least 6x6 minutes exposures with a random offset 
of 10 arc-seconds in between. Some details about the 50 
fields are given in Table 1. 

3. Data reduction 

The data were received by the end of November 1999 at 
lAP. The processing of the 50 fields was done using the 
TERAPIX^ data center facihties. The FORSl images are 
much smaller than the big CCD mosaics usually processed 
by TERAPIX, so we used a specific pipeline which in- 
cludes some IRAF scripts. Since the FORSl detector has 
four outputs, we simply consider each quadrant as an in- 
dividual device during all the preprocessing stage. 

The master bias was built using a standard median fil- 
ter. For the superfiat, since the CCD images do not show 
fringe patterns, the processing can be done in one step 
with a simple un-smoothed master-flat to correct for the 
relative pixel-to-pixel variations. However, the multiplic- 
ity of observations done in service mode provides images 
taken during different Moon periods and with different 
weather conditions which cannot be assembled to build 
up a single superflat. In our case, when all the images are 
used, the global master-fiat shows a strong and persistent 
gradient along the CCD (> 10%). It is therefore prefer- 
able to use daily master fiats by averaging only the images 
taken during the same night. The drawback is the small 
number of images which can be combined together. If it 
is too small, each individual frame has a strong relative 
weight which may results in correlated noise between in- 
dividual images and the master-flat. Therefore we decided 
to produce a special daily master-fiat for each image: each 
individual frame was flattened by using a master-flat done 
with the rest of the data taken during the night. 

The four quadrants are then put together and rescaled 
according to the gain of each output. If the sky back- 
grounds of the four quadrants are still different, we pro- 
ceed to another small rescaling accordingly This proce- 
dure turns out to work well in most cases. However, some 
images are not perfectly flatflelded and still show a small 
cross-shape discontinuity in the middle of the FORSl im- 
ages, between the borders of the quadrants. This area has 
been masked in the flnal co-added images when necessary. 

The stacking of each frame is done using standard 
IRAF procedure . The flnal error for alignment was less 
than one tenth of pixel. The flnal set of 50 stacked fields 



turns out to have a remarkable seeing distribution (see 
Figure 1). In fact, 90% of them are within our seeing spec- 
ifications, with an average value of 0.63". 



45 fields 



<seeing> = 0.631 



Fig. 1. Seeing distribution of the final VLT field sample used 
to measure the cosmic shear. 



.11 > 31.9^5 i 6.6:^ 



see http://terapix.iap.fr 
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Fig. 2. Histogram of the galaxy number density distribution 
in the VLT fields. The distribution shows a narrow peak at 
n = 35 gal arcmin"^. 



The photometric calibrations of the FORSl fields 
were done using the photometric standard star im- 
ages provided by the ESO staff at Paranal. The 
fields were chosen in a fairly broad sample of the 
Landoh's catalog (Landolt 1992): SA92, SA95, SAlOl, 
SA109, SAUO, MARK-A, PG0131-051, PG1313-086, 
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Table 1. List of the 50 VLT/FORSl fields. All have an exposure time of at least 36 minutes and cover a total area of 0.64 
deg^. The seeing given in this table correspond to the stacked image. 



Target Name 


RA (J2000) 


DEC (J2000) 


Seeing 


Ngal 


Iab Lim. 


vlt27 


00 


59 


28.1 


-00 


18 


28 


0.72" 


991 


24.5 


vlt28 


01 


31 


40.3 


-00 


22 


28 


0.54" 


1580 


24.9 


vlt29 


01 


59 


40.8 


-00 


03 


51 


0.49" 


1878 


25.1 


vlt30 


02 


28 


44.0 


-00 


03 


26 


0.54" 


1680 


24.9 


vltSl 


01 


00 


21.8 


-03 


15 


31 


0.57" 


1539 


25.0 


vlt33 


02 


00 


08.1 


-03 


00 


30 


0.44" 


1755 


24.9 


vlt35 


00 


59 


35.3 


-06 


10 


05 


0.73" 


1201 


24.8 


vlt36 


01 


28 


53.1 


-06 


01 


39 


0.68" 


1595 


24.8 


vlt37 


01 


57 


05.8 


-06 


05 


01 


0.90" 


1138 


24.5 


vlt39 


21 


30 


45.3 


-09 


58 


45 


0.76" 


1369 


24.8 


vlt40 


22 


04 


37.9 


-10 


15 


09 


0.71" 


1454 


24.9 


vlt42 


22 


29 


29.2 


-10 


12 


01 


0.72" 


1268 


24.7 


vlt43 


21 


30 


25.3 


-15 


11 


48 


0.55" 


1525 


25.0 


vlt44 


22 


02 


16.6 


-14 


53 


03 


0.64" 


1715 


24.9 


vlt45 


22 


30 


41.8 


-14 


54 


55 


0.46" 


1791 


24.8 


vlt46 


22 


01 


42.2 


-20 


10 


55 


0.65" 


1651 


25.0 


vlt47 


22 


29 


33.8 


-20 


14 


44 


0.51" 


1639 


25.0 


vlt48 


21 


30 


45.3 


-24 


53 


40 


0.63" 


1604 


24.8 


vlt49 


21 


58 


44.7 


-24 


57 


15 


0.62" 


1494 


24.6 


vlt50 


22 


30 


43.8 


-25 


01 


42 


0.55" 


1595 


25.0 


vltSl 


20 


59 


30.5 


-30 


18 


31 


0.62" 


1243 


24.5 


vlt52 


22 


00 


26.2 


-30 


01 


45 


0.65" 


1357 


24.8 


vlt53 


22 


31 


15.3 


-30 


07 


15 


0.55" 


1584 


24.7 


vlt54 


21 


29 


53.6 


-34 


51 


52 


0.57" 


1629 


25.1 


vlt55 


22 


00 


14.1 


-35 


30 


54 


0.53" 


1444 


24.6 


vlt56 


22 


30 


06.4 


-35 


10 


33 


0.83" 


1007 


24.5 


vlt57 


21 


28 


04.9 


-39 


49 


02 


0.55" 


1650 


25.0 


vlt58 


22 


00 


06.7 


-40 


04 


55 


0.49" 


1664 


24.8 


vlt59 


22 


29 


11.8 


-39 


36 


28 


0.70" 


1453 


24.7 


vlt60 


22 


59 


24.4 


-10 


01 


29 


0.47" 


1694 


25.0 


vitei 


22 


59 


24.2 


-15 


08 


47 


0.47" 


2115 


25.0 


vlt62 


22 


59 


01.8 


-19 


44 


03 


0.47" 


1769 


25.0 


vlt63 


22 


59 


39.5 


-24 


52 


51 


0.49" 


1502 


25.0 


vlt64 


22 


59 


56.1 


-30 


14 


27 


0.60" 


1285 


24.5 


vlt65 


23 


00 


44.3 


-34 


55 


26 


0.54" 


1675 


25.0 


vlt66 


23 


01 


24.8 


-40 


25 


20 


0.77" 


896 


23.6 


vlt75 


21 


28 


14.7 


-20 


07 


18 


0.56" 


1683 


24.9 


vlt76 


21 


32 


21.0 


-30 


25 


57 


0.63" 


1434 


24.6 


vlt77 


14 


59 


07.4 


00 07 54 


0.80" 


1414 


24.8 


vlt78 


14 


59 


03.2 


05 11 32 


0.50" 


1862 


25.0 


vlt79 


14 


59 


32.7 


10 13 19 


0.76" 


1142 


25.0 


vlt80 


15 


30 


17.5 


00 10 58 


0.60" 


1511 


24.8 


vlt81 


15 


29 


40.4 


04 54 10 


0.63" 


1405 


24.7 


vlt82 


15 


28 


59.7 


10 14 59 


0.59" 


1401 


24.6 


vlt83 


15 


59 


00.7 


-00 07 14 


0.87" 


1015 


24.3 


vlt84 


16 


03 


35.0 


05 10 46 


0.91" 


923 


23.8 


vlt85 


15 


56 


47.6 


10 17 28 


0.66" 


1275 


24.7 


vlt86 


16 


00 


30.1 


14 58 35 


0.78" 


1231 


24.6 


stis7ncw 


22 


28 


13.6 


-26 59 11 


0.76" 


860 


23.3 


stislOnew 


12 


28 


32.8 


02 10 05 


0.93" 


930 


24.0 
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PG1633+099, PG2213-006, PG2331+055. The photomet- 
ric zero-points were calculated rising magnitudes com- 
puted by SExtractor ^ and by rescaling the gains of each 
quadrant accordingly. We noticed a fluctuation of the zero- 
point of about ±0.05 magnitude along the observing pe- 
riod, probably produced by the important variation of the 
atmospheric conditions during the whole semester. 

The final photometric catalogs look sufficiently homo- 
geneous from field to field for our purpose. The magnitude 
histograms for each field have the same shapes and show 
a cutoff at Ij^B = 24.5 — 25.1. This cut-off defines our 
limiting magnitude. Five fields have a brighter magnitude 
cut-off at Iab = 23.5 — 24.5. These fields have the worse 
observational condition: seeing higher than 0.8" and high 
sky background due to the moonlight. Some fields have a 
higher number of co-added frames but it does not change 
significantly the final magnitude limit of the 50 fields. 

4. Measurement of the shear signal 

4-1- Catalogue generation and selection criteria 

The galaxies have been processed using our modified 
version of IMC AT software''. The modifications and the 
general analysis method have already been described in 
Section 3 of Van Waerbekc ct al. 2000. Here we remind 
shortly the selection procedure, paying attention only to 
the differences with the previous CFHT data analysis. 
The principal steps are the foUowings; 

— Masking the boundaries: the splitting procedure of ex- 
posures generates noisy thick strips surrounding the 
CCD, with lower signal-to- noise ratio than the rest of 
the image, and a grid-like pattern produced by the 
stacking procedure. We therefore mask automatically 
the boundaries of each image. 

— Masking the central cross-shape area: the four readout 
ports of the CCD delineate four well-defined quadrants 
with different readout noise, charge transfer efficiency 
and gain. At the junction of the quadrants, the sharp 
changes of CCD properties appear as a cross-shape 
area located in the middle of the CCD. In principle, 
this area is very small, but, like at the boundaries, 
the small shifts between each exposiire enlarge its size 
which thus appears as a broad region with smaller 
signal-to-noise ratio than the rest of the image. On 
the fields affected by this effect, the central cross-area 
is masked. 

— Masking objects: large areas around bright stars, sat- 
urated columns and diffraction patterns are masked 
on all images. They can produce spurious elongated 
structures which can be misinterpreted as shear signal. 

Likewise, we also mask areas around bright and/or ex- 
tended galaxies, elongated stripes produced by residu- 



als from asteroids, satellites or aircrafts. The fraction 

of objects removed by this masking procedure (bound- 
aries, cross-shape area and objects) varies from 10% to 
20%. 

Stars selection: we use a radius (r^) vs magnitude 
(mag.) plot to select unsaturated stars (we do not in- 
clude stars closer than one magnitude from the satu- 
ration level). The vertical box which encompasses the 
stars is carefully sized to avoid saturated stars and the 
crowed faint-magnitude areas of the rh — mag. diagram 
where galaxies and stars mix together. This domain 
contains very low signal-to-noise objects and is most 
of the time useless for the calibration. 
Fit of the stellar polarizability tensors: using these 
stars, we can map the PSF by fitting their raw shape to 
a third order polynomial. We then compute the quan- 
tities Pa = -p^ and -pTT^; 

Source selection: The starting sample of galaxies is 
extracted from the 50 VLT fields using SExtractor 
(Bertin & Arnouts 1996) and produces a raw cata- 
logue of 72500 galaxies. Objects at the edge of the 
fields and deblended objects are immediately removed. 
Galaxies with a half light radius smaller than the 
largest size of the star box are also removed. As pointed 
out by Van Waerbeke et al. 2000, close pairs of galax- 
ies may produce spurious elongated galaxies, conse- 
quently pairs with separation d < 10 pixels are also 
removed. We do not perform any magnitude selection, 
therefore the redshift distribution is broad and can be 
described by: 



p{z) = (3 



r((i + «)//?) f z 



— exp 



Zo 



(1) 



Using MAG_AUTO with a minimum radius of 3.5 pixels. 
® Made available by Nick Kaiser; 

http://www.ifa.hawaii.edu/~kaiser/ 



with a = 2, (3 = 1.5, = 0.8, which is consistent 
with a limiting magnitude Iab ~ 24.5 as given by 
Cohen et al. 1999. 
— Estimate of the shear: using the stellar fits, we compute 
the 'pre-seeing" shciar polarizability for all the selected 
sources as described in Van Waerbeke et al. 2000. 
Figure 3 shows the star ellipticities for all VLT fields 
before and after the PSF correction. One of the fields 
shows a seeing significantly alterated by strong wind 
and its image quality is too far from our specifications, 
this field (stis7new in Table 1) was removed. We ended 
up with 49 fields with in average 30 stars/field corre- 
sponding to a total number of 46941 galaxies, that is 
~ 21 sources/arcmin^. Note that this is less than the 
number given in Table 1 which was the number den- 
sity before objects selection through our shape mea- 
surement process. 

4-2. Shear measurement 

At this stage, the variance of the shear can be computed 
exactly as described in Van Waerbeke et al. 2000. How- 
ever here we have measured the signal in a slightly dif- 
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ferent manner. Instead of computing the variance of the 
shear by simply squaring the averaged shear per ceU, we 
directly removed the diagonal terms from the squared 
quantity, such that the computed variance is by defini- 
tion an unbiased estimate of the true shear variance (see 
Schneider et al. 1998). Therefore we do not need to sub- 
tract the shot noise contribution as in previous works, 
which was done using time-consuming monte-carlo ran- 
domizations. Let us call ea{6k) the ellipticity of a galaxy 
at a position 6k, and Wk its weight, calculated according to 
Van Waerbeke et al. 2000 (see Section 3.2 of that paper). 
An unbiased estimate of the variance of the shear 7^(0i) 




el el 



Fig. 3. PSF ellipticity before (left) and after (right) anisotropy 
correction. 

at location 6i can be directly obtained by measuring: 

2 N 

Eh\e,)] = "-^^^ (2) 

wkwi 

k^l 

The final product, the variance of the corrected ellip- 
ticity of galaxies as function of angular smoothing scales, 
is summarized on Figure 7. On this Figure we also plot 
the results of the previous detections. The angular scale 
corresponds to the diameter of a circular top-hat filter 
(we rescaled accordingly the angular scales of the previ- 
ous studies which were published for a squared top- hat). 
The 1-a error bars have been calculated from 100 real- 
izations of noisy catalogues by randomizing the position 
angle of each galaxy. 

The amplitude and the shape of the signal are in good 
agreement with former results. They are within the 1-a er- 
ror at all scales (note that the points are not independent). 
The agreement of these new measurements with previous 
ones provides a new independent indication that the signal 
is not produced by uncontrolled systematics. However, as 
for the CFHT analysis (see Van Waerbeke et al. 2000), we 
must analyze carefully the systematics in the VLT data. 

We checked the stability of this result with respect to 
the selection criteria by changing the star selection or the 



galaxy selection criteria. For the stars, we changed the 
box-size which encompasses the objects inside the vertical 
rfi — magnitude branch by increasing the magnitude range 
by 0.5 magnitude. For galaxies, we changed a few parame- 
ters, like the object definition (9 contiguous pixels having 
0.5 cr above the background or 6 contiguous pixels having 
1.5 cr above the background). We found that the ampli- 
tude of the variance fluctuates by ±7% on all scales but 
2.5 arc-minute for which the fluctuation is larger (±15%). 
These fluctuations are within the 1-a error bars, therefore 
selection criteria do not have a significant influence. 
Figure 7 shows an unusual behavior at 6 « 2.5 arc- 
minutes. The slope changes and the amplitude of the shear 
seems to increase and then to drop again on larger scales. 
This is a marginal perturbation within the error bars, but 
it turns out that it appears at the angular scale corre- 
sponding to one half the CCD size, so it could results 
from the masking procedure of the cross-shape area. We 
looked at this more carefully by computing the shear sig- 
nal inside vertical and horizontal strips of similar width 
to the mask. The results shown on Figure 4 (which is dis- 
cussed in the next Section) do not reveal any difference 
between the central regions {X = or Y = 0), where the 
two segments of the cross are located, and the rest of the 
image. We conclude that there is no evidence for system- 
atics generated by the 4-port readout configuration and 
that the increase of the signal is likely random error due 
to the lower signal-to-noise ratio on that scale. 

5. Analysis of systematics 

A critical issue regarding the measurement of cosmic shear 
signal is the understanding and handling of the various 
systematics. As in previous studies, we have taken special 
care of this point with this new VLT data set. 
Van Waerbeke et al. 2000 discussed various types of sys- 
tematics, attempted to provide some technical solutions 
to avoid part of them prior to use the KSB correction 
and provided a posteriori quality check on the shear sig- 
nal (see Sect. 5 of Van Waerbeke et al. 2000). We used 
similar controls for the VLT data: 

— CCD effects: The spurious signal produced by bad 
charge transfer efficiency, very bright stars, big galax- 
ies or asteroids are considerably reduced by the mask- 
ing procedure described in Section 4. Possible resid- 
uals from charge transfer efficiency can be estimated 
by looking at the variation of the ei and 62 compo- 
nents as function of the distance of objects with re- 
spect to the readout port, for each quadrant of the 
CCD. As shown in Figure 4, the corrected components 
(ei) and (62) do not show variations along the X or 
Y strips, except a small negative ei component ver- 
sus the Y strip for Y < 1000. We shall see in the 
following that it has a completely negligible contribu- 
tion to the variance of the shear. In fact, these plots 
do not show the significant systematic residual of the 
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Fig. 4. Average corrected galaxy ellipticity {ca) along verti- 
cal (Y) and horizontal (X) strips. The two components do not 
show any variation along the columns or the lines of the CCD. 
Moreover, the corrected (ei) component is well centered around 
zero and does not show the systematics negative value observed 
by Van Waerbeke et al. 2000 on CFHT data. The plots also in- 
clude the cross-shape area in the middle of the CCD where the 
four quadrants overlap {X = or Y = 0). We cannot see sig- 
nificant variations at those positions, which make us confident 
that the four port-readouts do not produce systematics. 



(ei) component which was observed on CFHT data 
(Van Waerbeke et al. 2000). Therefore it is hkely that 
the residual observed on CFHT data is intrinsic to the 
CFH12K camera. This also explains why it was not ob- 
served also by Kaiser et al. 2000 who only used UH8K 
data. 

residuals from the correction of the PSF anisotropy: 
FORSl has a remarkable image quality over the field. 
Because of the high quality of the optical design and of 
the small field of view of FORSl as compared to the 
UH8K and the CFH12K cameras, optical distortions 
are much smaller than in Van Waerbeke et al. 2000. 
However, active optics on the VLT could eventually 
produce a new and/or unexpected systematic effect 
not fully corrected with KSB. 



We tested the systematic residuals in exactly the same 
way as in Van Waerbeke et al. 2000. The results are 
summarized in Figures 5 and 6. The former shows that 
before the PSF anisotropy correction, the ellipticities 
of stars and galaxies arc correlated. In contrast, af- 
ter the correction the average galaxy ellipticity is zero, 
whatever the PSF anisotropy is. It shows that the cor- 
rection of the PSF anisotropy works very well even in 
the case of active optics and does not bias the corrected 
ellipticities of the galaxies. 

As pointed out in Van Waerbeke et al. 2000, the av- 
eraged ellipticity of galaxies binned, either with respect 
to the star ellipticity or to the CCD lines/columns as de- 
scribed above, is not a strong enough test of systemat- 
ics because it is still possible that the galaxy ellipiticities 
strongly fluctuate inside very small bins. It is therefore 
better to measure the variance (7^) instead of a simple 
average in bins. Moreover, the result found can be com- 
pared directly to the signal provided that, for each scale, 
the bin size is adapted to encompass a similar number of 
galaxies as in the top-hat filter used to measure the signal. 
If our corrections arc not contaminated by strong system- 
atics, at all scales this variance must be much smaller than 
the signal. 

Figure 6 shows the results. For each of the three panels, 
the filled circles with error bars show the VLT results as 
shown in Figure 7, and the dashed lines show the 1-a er- 
ror bars obtained from 100 randomizations of the position 
angle of galaxies. The open circles in the top panel show 
(7^) measured in bins of galaxies sorted according to the 
strength of the anisotropy of stars (cither e* or corre- 
sponding to the two sets of open circles). The open circles 
in the middle panel show (7^) with the galaxies sorted 
according to the X and Y positions on the CCDs, (again, 
each set of open circles correspond to galaxies sorted either 
according X or to Y). It is clear from this plot that the 
slight negative component (ei) observed for Y < 1000 in 
Figure 4 is not strong enough to produce a significant sys- 
tematic compared to the signal. These two panels demon- 
strate that the residuals are always within the l-cr error 
and fluctuate around zero. We conclude that these resid- 
uals are not responsible for the signal we detect. 
It is interesting to estimate how strong the systematic 
effects could be if we neglected to correct the shape of 
the galaxies for the PSF anisotropics. This is shown in 
the bottom panel of Figure 6 where (7^) is computed in 
the same way as for the top panel, but with galaxies un- 
corrected for the PSF anisotropy. The uncorrected signal 
shows a significant offset which is almost insensitive to the 
angular scale (because the PSF is roughly constant over 
each CCD). Therefore, the shape and the amplitude of the 
signal are unlikely to be produced by residuals from the 
PSF anisotropy. Moreover, we see that a complete lack 
of anisotropy correction produces a systematic of similar 
amplitude as the real signal (and not much higher). There- 
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Fig. 5. Average galaxy ellipticity (ea) compared to the average 
star ellipticity ^e^*'"'^. The two plots at the top show the aver- 
age galaxy ellipticity before correction for the PSF anisotropy. 
There is a tight correlation between the galaxy and star el- 
lipticity, as expected. The two plots at the bottom show the 
same components after the correction: galaxy and star elliptic- 
ities are no longer correlated and the two components of the 
average galaxy ellipticity are centered on zero. This demon- 
strates that the correction works well and does not produce 
any bias. 



6. Discussion 

The VLT data confirms the detection of a significant weak 
distortion signal on angular scales between 0.5 to 5 ar- 
cminutes. Its amplitude and its shape are similar to those 
announced previously on these angular scales by four inde- 
pendent teams. The study of systematics does not reveal 
any bias and reproduces similar trends as for the CFHT 
data but on a much larger sample of uncorrelated fields 
and on a more homogeneous sample (only I band, narrow 
seeing distribution and depth). 
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Fig. 6. Analysis of possible systematics on the VLT fields com- 
pared to the signal. Filled circles with error bars correspond to 
the cosmic shear signal; the dashed lines shows the 1-a error 
bars obtained from the 100 randomizations of VLT galaxy ori- 
entations. The top panel shows (7^) (open circles) measured 
on the galaxies sorted according to the PSF of stars (see the 
two bottom panels in Figure 5). (l'^) is measured in bins of 
galaxies, which has a number of galaxies comparable to the 
number of galaxies present in a top-hat filter of radius 9. The 
middle panel shows (7^) (open circles) measured on galaxies 
sorted according to their Y and Y location on the CCDs (see 
Figure 4). Finally, the bottom panel shows (7^) (open circles) 
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It is interesting to investigate what constraints on cos- 
mological models wc could provide from the cosmic shear 
surveys completed so far 

by Van Waerbeke et al. 2000 (CFHT), Bacon et al. 2000a 
(WHT), Wittman ct al. 2000 (CTIO), Kaiser et al. 2000 
(CFHT) and by adding the VLT data (this work). The five 
data set have been observed either in R or in I-bands, at 
roughly the same depth (/ « 24.0), so we can assume that 
the average redshifts of the sources are almost the same 
and should be close to the value inferred from the deep red- 
shift survey carried by Cohen et al. 1999: zq « 0.8 with 
the broad distribution given by Eq.(l). Figure 7 shows 
some current cosmological models with the present-day 
data. The non-linear evolution of the power spectrum of 
density fluctuations has been taken into account following 
the prescription given by Peacock & Dodds 1996^. One 
can see that the cluster normalised models fit very well 
the observations, at least on scales ranging from 0.5 to 10 
arc-minutes. However, this plot does not really illustrate 
the constraints on both Qq and erg we can expect from 
measurements of the variance of the shear. A more reli- 
able study consists in exploring a very large set of models 
in a (riojCg) space. As we can see from Figure 7, and as 
noted before (Bernardeau et al. 1997), the dependence on 
the cosmological constant of the variance of the shear is 
rather weak, and it is not worth to include A as a free 
parameter in this analysis. 

Let us consider all the five cosmic shear results simul- 
taneously. Since they provide independent samples, we 
can use one single measurement point for each of them 
and perform a simple minimization in the (fiefs) 
plane. From each cosmic shear measurements we choose 
the point which has the best signal-to-noise, and we 
avoid the large scale measures in Kaiser et al. 2000 and 
Wittman et al. 2000, as they are likely affected by fi- 
nite size effects which tend to increase the error bars 
(see Szapudi &: Colombi 1996 for a general discussion, or 
Bernardeau et al. 1997 for a specific application to cosmic 
shear surveys). We extracted five triplets containing the 
scale, the variance and the 1-cr error, {0i, j'^{9i),S^'^{9i)), 
out of the literature (the results reported on Figure 7) and 
computed: 



3-a confidence level contours. The best fitted models can 
be described by the empirical law: 



E 

i=l 



(3) 



where {'y^)ei is the predicted variance for a given cosmo- 
logical model. Wc computed it for 150 models inside the 
box < < 1 and 0.2 < as < 1.4, with F = 0.21, 
A = 0, and zo = 0.8. With 5 data points and two free 

parameters, the has 3 degrees of freedom. The result is 
given in Figure 8. The grey scales indicates the 1, 2 and 



® Following Peacock's advice, we do not use anymore the co- 
efficients given in his textbook, "Cosmological Physics" which 
turn out to be less accurate than in their paper. 



as ~ 0.59+^; 



-1-0.03 Q-0A7 



03 "0 



(4) 



in the range 0.5 < 6 < 5 arc-minutes. This is in re- 
markable agreement with Jain & Scljak 1997 who pre- 
dicted non-linear scales, and this is very 
close to the cluster normalization constraints given in 
Pierpaoli et al. 2000 (for closed models and F = 0.23): 



as ^ 0.495t°;0r, flo"' 



60 



(5) 



This law overplotted on Figure 8 shows a remarkable 
agreement between these two approaches. 

Our analysis is still preliminary. We have only five in- 
dependent data points spread over a rather small angular 
scale. We also assumed the peak in the rcdshift distribu- 
tion of the sources to be zq = 0.8 with the broad distri- 
bution given in Eq.(l). This is a reasonable assumption 
on the basis of the spectroscopic survey carried out by 
Cohen et al. 1999, but it is still uncertain and needs fur- 
ther confirmations. We also neglected the cosmic variance 
in the error budget of the cosmic shear sample. Although 
it does not affect seriously the VLT data which contain 50 
uncorrelated fields and the Bacon et al. 2000a data point 
(because they estimated the cosmic variance using a Gaus- 
sian field hypothesis), the three other measures are prob- 
ably affected. However, numerical simulations already in- 
dicate that cosmic variance should only increase our error 
bars by less than a factor of two, according to the survey 
size (see Van Waerbeke et al. 2000). We expect to have 
much better constraints on the redshift and the clustering 
of sources once the VIRMOS redshift survey will be com- 
pleted (Le Fevre et al. 2000). On the other hand, a mea- 
surement of the skewness of the convergence will break the 
degeneracy between Qq and as (Bernardeau et al. 1997, 
Van Waerbeke et al. 1999). 

7. Conclusion 

Wc have confirmed the cosmic shear signal detected on 
scales ranging from 0.5 to 5 arc- minutes using for the first 
time a large sample (50) of uncorrelated fields obtained 
with the VLT UTl/ANTU. The service mode available 
on this telescope gave an unprecedented high quality and 
homogeneous data set. The fields are spread over more 
than 1000 square-degrees, which minimizes the noise pro- 
duced by cosmic variance. The amplitude and the shape 
of the shear are similar to those measured from other tele- 
scopes which permits to make a strong statement about 
the robustness of the signal with respect to various sources 
of systematics. 

Assuming the signal is purely cosmic shear (that is we 
neglect the possible intrinsic shape correlation), we used 
the four other studies published so far to infer first con- 
straints on cosmological models. The cosmic shear surveys 
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Fig. 7. (7'^) as function of the angular scale {0 is the angular 
diameter of a circular top-hat). The filled squares are the VLT 
data presented in this work. We also plot the previous detec- 
tion: Van Waerbeke et al. 2000 (vWME-|-), Bacon et al. 2000a 
(BRE), Wittman et al. 2000 (WTK) and Kaiser et al. 2000 
(KWL). These measurements are compared to some current 
cosmological models with galaxies following the redshift distri- 
bution Eq.(l) with zo — 0.8. For all the models we choose a 
CDM power spectrum with F — 0.21, and f2o = 1, A = 0, 
o-g = 1 (short dash); Qq = 0.3, A = 0, erg = 1.02 (dot 
long-dash); Slo ~ 1-0, A = 0, erg = 0.6 (dot) and Q,o — 0.3, 
A = 0.7, erg = 1.02 (solid). The models have been computed 
using the non-linear evolution of the power spectrum given by 
Peacock & Dodds 1996. 



provide constraints which are in remarkable agreement 
with those from cluster abundance analysis. As compare 
to it, the variance of the shear is a direct measure of the 
combined parameters flo,as- As soon as the residual bi- 
ases are well controlled, and the redsfhit of the sources 
known, the measurement of (7^) will naturally converge 
to the exact D,q, erg value with an accuracy that will only 
depends on the accumulation of measurements on uncorre- 
lated fields of view. This study shows the great potential of 
cosmic shear for cosmology and what we can expect from 
future wide fields surveys. In particular, the skewness of 
the convergence, which is insensitive to cts, will appear as 
a vertical constraint on Figure 8 therefore breaking the 
flo-(Ts degeneracy. 
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